GRBs: when do blackbody spectra look like non-thermal ones? 
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ABSTRACT 

We argue that a nonthermally looking spectrum of a gamma-ray burst (GRB) can 
,^ • be formed as a superposition of a set of thermal blackbody spectra. This superposition 

[t . may be done by time integration which is present even in 'time resolved' GRB 

1^^ i spectroscopy. A nonthermal spectrum can be obtained also by the space integration 

CN i which should take place unless all the emission comes from a plane front moving 



0^ 



(N 



exactly towards the observer. We propose a model of the gamma-ray burst spectrum 

^ ' formation based on this idea. This model allows the GRB radiation to be optically 

l>^___ ■ thick and to have higher values of baryon load. Thus the latter is limited by the energy 

cn ' considerations only, and not by the condition of a small optical depth. 

(N 

O 

o\ 

^* ' Subject headings: gamma-rays: bursts — gamma-rays: theory — radiation 

T^ • mechanisms: thermal 

o 

•^ ■ 1. Motivation 

>'■ 

K> , Gamma-ray bursts (GRBs) still remain an unresolved mystery of modern astrophysics in spite 

H I of recent progress in the observations of their X-ray, optical and radio counterparts. Not only the 

nature of internal engine, but even the mechanism of gamma-ray emission is unclear. Studying 
the spectra of GRBs is one of the keys that can unlock this great mystery in future. 

Observations of the GRB spectra (Band et al. 1993) show that, in general, they are well 
described by a low-energy power law with the exponent a, being exponentially cut of at -B ~ Eq, and 
by a high-energy power law with the exponent (3. Though the values of (a, /3, £"0) can be different 
for individual bursts, they usually are in the range of (~ —1.5 . . . 0.5, —3 ... 2, 100 . . . 200keV). 

Note that in this paper we consider the photon spectrum N{E) or N{v), the differential 
energy flux density Fy = hvN{hu), and vFy distribution. By default, all the power indices in this 
paper refer to N{E). 
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The power-law appearance of the spectra can possibly be explained by the hypothesis of their 
nonthermal origin. The synchrotron shock mechanism (Tavani 1996), where the GRB emission is 
produced by an optically thin relativistic plasma in a weak magnetic field, is one of those models 
which give a good agreement with observed spectra. Cohen et al (1997) find that the low-energy 
spectral index a in the time-integrated of GRB is usually in the range from —2/3 to —3/2 as 
predicted by the synchrotron shock model. The limits of this range correspond to the synchrotron 
spectra of instantaneous sample of electrons and the one integrated over their radiative decay 
(Rybicki k Lightman 1979). 

However, Crider, Liang & Preece (1997) have shown, on the basis of the analysis of the 
time-resolved spectra of 99 GRBs, that neither the synchrotron shock nor the simple inverse 
Compton mechanism can explain the instantaneous GRB spectra and their evolution: the 
time-resolved spectral slope a is often outside the limits of the synchrotron model and does not 
change monotonically with time, as the inverse Compton model predicts. 

While these models of gamma-ray bursts (which generally fit the observations) have some 
difficulties to match them in detail, we can present here a blackbody model that should be at least 
not worse than the other current ones. 

The conflict of the optically thick model for GRBs with observations was discussed already 
by Paczyhski (1986) and Goodman (1986). Paczyhski (1986) mentioned: 'The observed spectra 
are averaged over large fractions of a second, and this may be responsible for the shallow slope of 
the low energy part of the spectrum'. The problem was raised recently by Band &: Ford (1997). 
They have posed a question 'whether burst spectra are narrowband on short time-scales'. So, the 
question is: are the observed broadband GRB spectra formed by time integrations of an evolving 
quasi-blackbody instantaneous spectrum, or not. Band & Ford (1997) found no evidence for 
narrowband emission down to 1 ms time-scale. In the present paper we consider time-scales that 
are shorter for an observer. 

It is well known, that assuming high values of Lorentz factor F of the GRB ejecta is necessary 
to solve the compactness problem (Guilbert, Fabian & Rees 1983, Paczyhski 1986, Goodman 1986, 
Krolik & Pier 1991, Rees & Meszaros 1992, Piran 1996). The typical time-scale of the variability 
of the gamma-ray emission At ~ 10^^ seconds implies the size of the emitting region R < cAt, 
as small as ~ 10^ km. The enormous number of gamma photons in such a small volume should 
produce electron-positron pairs which make the emitting region optically thick. This conflicts 
with the observed nonthermal spectra unless one supposes that the emitting region moves towards 
the observer at a relativistic speed with Lorentz factor F, then its size would be F^cAt, and the 
optical depth correspondingly smaller. 

We propose an important supplement to this solution of the compactness problem. In our 
version, the relativistic motion is still required, in order to provide the formation of an integrated 
spectra from an ensemble of the thermal ones. 

It is known that a sum of different thermal blackbody spectra can produce a power-like 



spectrum looking as a nonthermal one. It happens, e.g., in the classical case of the Shakura- 
Sunyaev thin accretion disk (Shakura & Sunyaev 1973). As shown in this paper, a similar approach 
can provide an analogous result in the case of a relativistically moving emitter. 

Evidently, in any realistic situation the spectrum produced by an optically thick body is 
never a pure blackbody, because of opacity (and hence emission) dependence on wavelength, the 
effects of sphericity (see Mihalas 1978) etc. For us the black body is just a 'toy' model which is 
however far enough from the spectra of an optically thin plasma, invoked by others for explaining 
GRBs. Ryde & Svensson (1999) consider another basic model (a non-thermal one) and show that 
the observed spectra result from the time integration. Our approach is more radical than that. 

By the spectrum formation model presented in this paper we do not introduce a new physical 
model of gamma-ray bursts. We simply point out the fact that the observed non-thermal spectrum 
can be produced by an optically thick body. The assumptions needed for this seem not to be 
very unnatural. If such a picture can be worked out as a physical one (not only the 'toy' model), 
then new classes of GRB models become possible, producing 'dirty' fireballs, e.g. by the neutrino 
annihilation (Goodman, Dar, & Nussinov 1987). On the GRB models with a moderately high 
baryon load see Woosley (1993), Ruffert et al. (1997), Fuller & Shi (1998), Fryer & Woosley 
(1998), Popham, Woosley, & Fryer (1998). 



2. The model of spectrum formation 

Let us assume that the emitting surface is moving towards the observer with F ~ 10^ - it 
can be an expanding shell, or a blob, or a 'bullet', or an 'internal shock' (e.g. Piran 1998) - and 
producing at each instant a pure blackbody spectrum (which has a resemblance to the real spectra 
of optically thick plasmas). 

Due to the well known effect, if the emitter is moving towards the observer with the velocity 

V corresponding to F = (1 — t;^/c^)~^'^ then the emitter and observer time-scales differ by a factor 
of 2F2 (e.g. Rees & Meszaros 1992, Shaviv & Dar 1995, Piran 1998, Dar 1998). Here and below 
we assume that all clocks are synchronized in the observer's rest frame, i.e. the effect is purely 
kinematical (see Fig.|^), moreover it is Galilean, not truly relativistic (in the sense that Relativity 
plays no role in this effect). The Lorentz factor F is here simply a measure of the deviation of 

V from c, and nothing else. The difference of the emitter and observer time-scales means that, 
for example, r = 10 ms, the time of integration by an observer, corresponds to r' ~ 5 hours of 
emission time (Fig.|l|). During this long time the emitting object can expand and cool significantly, 
so the spectra it produces in the beginning and at the end of the observation interval r can differ 
drastically. Therefore, the observed spectrum is formed by an integration of some cooling sample 
of instantaneous spectra. 

For simplicity, we assume that the temperature T and the area A of the emitting object 
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Fig. 1. — The space-time diagram for the emission of a shell expanding with velocity /3c. The 
observer detects the duration r for the pulse emitted by the shell during the interval t'^ . . .t'^. 



change with time as described by the following power laws: 

T = To{t/to)-' = To{t'/Q~'; A = A^it/t^Y = A^{t' /t'^Y . (1) 

Here and below the primed time variables will refer to the emission time, while the non-primed 
ones - to the detection time. 



2.1. Analytic treatment 

Let us consider a set of arbitrary elementary spectra. If the members of the set have a 
parameter distributed according to a power law, then the integration of elementary spectra leads 
quite often to the formation of a power spectrum. Let us show this with a simple example (Fig. |2|). 

We would like to denote the elementary (instant) spectrum as n(E), and the resulting 
(integral) one as N(E). The spectra will be integrated in time from to to ti = to + ''"> where 
ti > to. 

1. Let the elementary spectrum be n{E) oc E~^ (Fig. ^) in the high energy part of the 
spectrum {E > £"0) and constant if {E < Eq), where Eq evolves in time like t~^ and at ti 
reaches the value Ei = £'o(*i/*o)^^- Then the observed integral spectrum should be 

ti r E~^, E>Eq 

N{E)= j A{t)n{E,t)dt^\ E-l^-^l\ Ei < E < Eq , (2) 

to [ const, E < El 

i.e. have two power-law parts the harder of which reflect the high-energy tail of the 
elementary spectrum and the softer accounts for the elementary spectrum evolution. 

2. Let us now examine the case of a stepwise elementary spectrum (Fig. 0b) described by a 
Heaviside 0-function: n{E) = Q{Eq — E), where £"0 evolves as in the previous example. 
Then the integral spectrum should be 

r 0, E> Eo 

N{E) = I r^E-^, Ei<E<Eo , (3) 

[ const, E < El 

i.e the Heaviside step is smoothed into a power function. The discontinuity of the above 
N{E) at the point Eq is an artifact of the approximation: in fact, there is not an exact 
power function, but very close to it if ti S> to as supposed. 

3. So we have shown that the integration of both power and stepwise spectra leads to a 
power-law behaviour between Ei and Eq. The elementary spectrum with Plank (Wien) 
high-energy tail (Fig. ^) lies between the power and stepwise cases, it is not as steep as 
the stepwise one but steeper than the power one. So it is natural to expect a similar result 
(power-law behaviour) for the integral spectrum. 
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Fig. 2. — The three cases of simple elementary spectra: with the power (a), stepwise (b), and Plank 
(c) high-energy tails. 
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2.2. The integration of the blackbody elementary spectrum 

Now we have come to the time integration of the blackbody Plank spectrum. 

"'^■" = -^"> exp[E/TW|-l - <"> 

Here we measure E and T in the same units, say, Tq, and let us measure time in units of t^, so 
instead of 

Ait) = A„(l-y. n*) = r„(l)" (5) 

we have simply 

A{t) = ^ot", T{t) = t-'. (6) 

The observed integral spectrum is 

N{E) = / dtA{t) = Ao / dt . (7) 

Ji exp[E/T{t)\ - 1 Ji exp{Er) - 1 

Introducing y = Et , we rewrite this as 

N{E) = Ao / dy ' . (8) 

9 Je exp(y) - 1 

From this general expression we derive the asymptotic cases. 

1. The most interesting case is when E < 1, that is i;^ < fcTg in standard units, and Eti ^ 1. 
One should remember that ti is always greater than unity, so the latter inequality is true 
when E is not too small. Then we find, replacing the lower integration limit by zero, and 
the upper one by infinity, 

N{E)c,Ao / dy ^ (9) 

Jo exp(|/) - 1 

The value of the integral is not interesting for us now. Thus, we produce a power-law 
spectrum with the exponent 2 — (a + l)/9. Say, for a = 2 and 9 = 3/4 we find the spectrum 
N{E) ~ z^~2. For cr = 2 and 6^ = 1 we find the spectrum N{E) ~ E''^ (flat F^ ~ z/O), etc. 
See the numerical examples below. 

2. When E <^ 1 and Eti <^ 1, we have y <^ 1, so we are in the Rayleigh-Jeans (RJ) regime, 
exp(y) — 1 ~ 2/) s-iid 

N{E) ~ Ao- / dyy^'^+^y^-^ oc E . (10) 



3. For high frequencies, E > 1 and Eti >> E > 1, the flux reduces to 

Je o 

(11) 
So, here, in the Wien regime, for any a and 9 we have in standard units 
N{E) ~ E^eyi^{-E/kTQ) [contrary to Nb{E) ~ E"^ e^-^i-E/kTo) for the blackbody 
of temperature Tq] . 

2.3. Numerical Examples 

Fig. |3| presents the results of numerical integration of the 4 cases of elementary spectra with 
various model parameters. It also illustrates the correctness of the above analytical estimates. 
One can compare these spectra with fig. 4 from Chiang & Dermer (1998) where a similar time 
integration is done, but in a different physical situation. 

For a fixed pair of a and 9 the spectrum consists of two power laws and one exponential 
(Wien) high-energy part. Therefore, it is in some sense similar to the Band (1993) function which 
also has two power law parts, so it can be expected to fit the observations as well. 

The moderately high energy part [Ei < E < Eq) has the power law spectrum with the the 
exponent depending both on the cooling (9) and expansion (a) laws. The dynamical range, i.e. 
the spectral width of this part is Ei/Eq = (ti/io)~^ ^ 1- Of course it depends on the integration 
time T = ti — to ~ ^1) ^-ncl should be smaller for high temporal resolution. This can be a serious 
test of the present model. 

The highest energy part [E > Eq) represents the exponential breakdown, which may be 
observed or not, depending on the value of Eq. Also, for such energies, there may exists some 
other (optically thin?) radiation mechanisms which can provide more intense emission than the 
proposed blackbody one. 

The low energy [E < Ei) part of the spectrum in our model should have only one possible 
value of the slope. This is clear from the analytical considerations: a = 1 (see eqn. |lO|). This a 
is close to be consistent with many observations (Crider et al. 1997), but the observed variety of 
spectra is much richer, than the simple RJ case, and there are claims that some GRB's do show 
here the spectra predicted by synchrotron model (Cohen et al. 1997). We can demonstrate, that 
with a small sophistication our blackbody model can reproduce those spectra as well. 

Introducing new parameters is usually a means to improve a fit, but also makes the latter 
physically less reliable. In what follows, we will keep one parameter constant, let us take o" = 2, as 
the most natural choice. Instead, we can introduce a physically motivated additional parameter 
/hard ^s the fraction of the time when the value of 9 is constant, assuming that after some time, 
/hard''") the temperature power law changes. In the examples below, for the fixed a at constant 
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Fig. 3. — The integrated spectra: a) step function with Tq = lO^eV, a = 2, i 
b) blackbody with Tq = 10^ 10^, lO^eV, cr = 2, Q = 0.75; c) the same for Tq 
0- = 2, 1.25, 0.5, e = 0.75; d) the same for Tq = lO^eV, a = 2, 6* = 1.5, 1, 0.75 
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value 2, we allow the value of 9 to change a bit. For illustration, we have taken GRB spectra 
given in (Cohen et al. 1997) in the form of postscript files and superimposed them onto our 
fits. The spectra from (Cohen et al. 1997) are all integrated in time, it would be better to have 
time-resolved spectra. But in any case it is not possible to have time resolution better then 1 ms 
and for the illustration of our idea the spectra used are quite good. As shown by Figs. |^, |5|, ^, 0, 
our black body model can provide good fits for the GRB spectra which were claimed to give 
evidence for synchrotron radiation. 



3. Discussion and Conclusions 

We have assumed that at each moment the spectrum of the gamma-ray burst emission is 
close to the black body one. After the integration in time over the typical temporal resolution of 
the observations it produces a spectrum which can be similar to the observed 'non-thermal' GRB 
spectrum. In reality, both the instantaneous spectrum and its true time evolution can deviate 
significantly from our simplified assumption. So in reality one can have a much richer variety 
of observed spectra. In our work, we wish only to point out a simple fact: that the observed 
non-thermal spectrum can be produced by an optically thick expanding body under fairly natural 
assumptions 

We have in mind the following picture. The central engine of GRB operates on a space scale 
like 10^ cm (the size of a neutron star or a stellar mass black hole). It produces shells or bullets of 
matter moving with the speed which is only one millionth slower than the light speed c, i.e. we 
assume F ~ 10^. The high value of F is needed in any case for cosmological GRBs in order to 
solve the compactness problem (see e.g. Piran 1998 for all refs). But the standard picture invokes 
the high F in order to make the fireball transparent only if it is clean (without baryon load) : they 
have the optical depth r^^ going down ~ F'^+^, if the /3 ~ 2 is the index of the power spectrum at 
the hard tail A^(^) oc E'l^. 

So in standard picture, and in our picture as well, if we see a pulse of GRB lasting 1 ms, the 
size of the shell should have grown from 10^ cm up to 10^ light milliseconds = 10^ light seconds 
~ 10^"^ cm, since R ~ 2F^c x 1 ms. Now one can only start speculating, where do the next pulses 
of GRB come from. These can be internal or external shocks (Piran 1998), or light reflections 
(Shaviv & Dar 1995, Drozdova & Panchenko 1997), etc. However, there are arguments (Fenimore, 
Ramirez, & Sumner 1997) that one shell expanding forever is not able to produce GRB pulses 
which only show a slight 'hard to soft' evolution for hundreds of pulses. Already for the first pulse 
the shell had to expand from ~ 10^ cm to -^ 10^^ cm. Therefore, a model where a central engine 
repeats shooting shells or bullets for the whole duration of the GRB is preferred (see also Dar 
1998). 

Thus, we have R like 10^ cm and to ~ 3 • 10~^ sec, and ti (or r) like 10^ second, so 7 orders of 
magnitude for the dynamical range of a power-law spectrum in our model is quite plausible. This 
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Fig. 4.— a) GRB920406: = 0.75, 9s = 1.0, IgTo = 5.45, h/to = 1.2 • 10^, /h^rd = 0.0115; b) 
GRB920622: 6 = 0.75, Os = 0.955, IgTo = 5.62, ti/to = 1.7 • 10^, /hard = 0.05 



12 




100 

Energy (KeV) 




100 

Energy (KeV) 



Fig. 5.— a) GRB92n23: 9 = 0.75, 6s = 1.125, IgTo = 5.28, h/to = 8.7- 10^, /hard = 0.01 b) 
GRB930201: 9 = 0.75, 9^ = 0.975, IgTo = 5.47, h/to = 1.3 • 10^, /hard = 0.014 
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Fig. 6.— a) GRB930506: 9 = 0.74, 6s = 0.967, IgTo = 6.0, h/to = 1.1 • 10^, /hard = 0.014 b) 
GRB940206: 9 = 0.74, 9^ = 0.975, IgTo = 5.65, h/to = 1.1 • 10^, /hard = 0.014 
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covers the range from keV to GeV. The evidence for hard, TeV, emission associated with GRBs 
remains inconclusive (Padilla et al. 1998). Only if the extremely hard TeV photons are detected, 
as suggested by Totani (1998), then one should invoke a truly non-thermal emission mechanism. 
We do not say that our shell has the thickness in the end like ~ 10^^ cm, it must be geometrically 
thin, so more dense - it is optically thick. But its radius R is of course ~ 10^^ cm, with clR -C R. 
It can be loaded with baryons to some extend, not violating the energy limits of course (Krolik & 
Pier 1991). This is good if one has something like stripping the surface layers of neutron stars 
(Blinnikov et al. 1984; Eichler et al. 1989; Ruffert et al. 1997). Reaching the size ~ 10^^ cm our 
shell (or bullet) has expanded and cooled enough to become transparent in the end. The shell 
traveled this distance 10^ seconds according to our clocks, but one should not forget that it kept 
running almost with speed of light, the light that it had produced. So the difference in time of the 
beginning of the flash, that we see on Earth, and its end is only 1 millisecond. While our shell is 
still very near the centre, the engine has shot already the second shell (or bullet, Fig. ^, then the 
3rd one, ..., the 100th, etc. If the total GRB duration observed on Earth was a few seconds, all 
the shots of the central engine were done when our first shell was like (few seconds/ 10^ seconds) 
smaller than in the end, so its radius was like ~ 10^^ to ~ 10^^ cm. At this time it was very 
optically thick, but one should not forget that it moves so fast, that the light of the 2nd, 3rd, ..., 
100th etc. shells can reach the first shell only after the first one is far away, ~ 10^*^ cm from the 
centre, and absolutely transparent. 

If instead of shells we have bullets, moving at some small angles to us there is no problem of 
transparency. They can cool down and become small solid bodies (this is perhaps not probable, 
since they must be heated up by ISM). 

In reality, not only time, but also space integration takes place. As shown by Rees (1966), 
(see also Drozdova &: Panchenko 1997, Sari 1998) in the case of an expanding emitting shell 
an observer simultaneously detects radiation produced in different moments of time (thus, with 
different temperatures) on the ellipsoidal or egg-like surface. The integration over this surface can 
give the same effect as the integration over time done in this paper, but we do not perform this 
here because the result strongly depends on the unknown geometry of the emitting surface. 

To conclude, we found that a variety of observed 'nonthermal' GRB spectra can be well 
reproduced by the time-integrated emission of a black-body spectra. The most critical test of our 
model can be the discovery of the temporal resolution dependence of the power spectrum range 
(here Ei . . . Eq). However, it can be smoothed by a space integration. The main advantage of the 
proposed model is that it allows the baryon load to be limited not by the optical thickness, but 
by energy considerations only (one cannot accelerate too much baryons because of their high rest 
mass) . 
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Fig. 7.— GRB940302: 6 = 0.75, ds = 0.96, IgTo = 5.45, ti/to = 7.8 • 10^, /hard = 0.012 
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Fig. 8. — The space-time diagram for the emission of two shells. Though the second shell is 
shot when the first one remains optically thick, the observer sees the radiation of these two shells 
separately because of their ultra-relativistic velocities. 
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